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Do dwarf galaxies form in tidal tails? 
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ABSTRACT 

The formation of tidal dwarf galaxies (TDG) inside tidal arms of interacting disk 
galaxies has been studied with iV-body and iV-body/SPH simulations at different res- 
olutions. In pure A^-body simulations no bound objects are formed at high resolution. 
At low resolution bound objects can form in tidal tails in agreement with previous 
work. We conclude that tidal dwarf galaxies are not likely to form by pure collisionlcss 
collapse in tidal tails. However, the presence of a sufficiently massive and extended 
gas component in the progenitor disk supports the formation of bound stellar objects 
in the tidal arms. Our results clearly favor a dissipation supported scenario in which 
the formation of TDGs is induced by the local collapse of gas which then triggers the 
collapse of the stellar component. 
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1 INTRODUCTION 

1.1 Observations of Tidal Tails in Interacting 
Galaxies 



Interacting galaxies are very complex and highly dynamic 
systems. They can be comprised of spiral galaxies as well 
as elliptical galaxies which possibly merge and change mor- 
phological type. If the mass ratio of spiral mergers is in the 
range 1 : 1 to 4 : 1, from theoretical as well as observational 
point of vie w the systems are l ikely to e volve into ellipti- 
cal g alaxies (jGenzel et al.l l200ll . but see iRobertson et al.l 
120061 for a counter example) . For those ma ss ratios the rem- 
nant properties resemble elli ptical galaxies dNaab fc Burkertl 
120031 : iNaab fc Truiilloll2006l ). Typically, mergers take place 
in field and group environments where the velocity disper- 
sion is of the order of the escape speed of the galaxies. How- 
ever, me rging has also been observed in probably u nrelaxed 
clusters l|van Dokkum et al.lll999l ; iTran et al.ll2005l ). 

Mergers between spheroidal galaxies are not spectacu- 
lar. As these systems are mostly gas poor and dynamically 
hot the merger is neithe r accompanied by a starburst nor by 
distinct tidal features (Ivan Dokkum et al 1 ll999l ; iBell et al.l 
l2006l : lNaab et all 120061} . Hence the formation of new dwarf 
galaxies inside tidal features is very unlikely. 

If spiral galaxies with gas and dynamically cold 
discs are involved in a merger the systems show strik- 
ing features like starbursts, AGN activity jGenzel et alj 
1 19981 : ISanders fc Mirabellll996l ; iDasvra et al.ll2006l ) and long 
arms of stellar and gaseous material which extend from 
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the parent galaxy to sometimes hundreds of kiloparsecs 
away. Some famous examples are NGC 4676 ('The Mice'), 
NGC 4038/39 ('The Anten nae'), IRAS 19254-7245 ('The 
Superantennae') (s e e e.g. iHibbard fc van Gorkoml Il996l: 
IWilson et all l2000l: IHibbard et all l200ll: iFabbiano et al. 



1 2003) : IHibbard et alj|2005l; iMirabel et al1ll99ll ; IVanzi et al 
2002; ICharmandaris et al l l2002i . and references therein). 
These arms are produced during the close encounters of the 
galaxies when strong tidal for ces expel the material fro m the 
discs of their parent galaxies l|Toomre fc Toomrell 19721 ). Al- 
ready HidckJ HHl) 

suggested that inside these tidal tails 
ne w stellar systems could be forming. This idea was revived 
bv lSchweizerl l| 19781) who found a low surface brightness ob- 
ject slightly beyond the tip of the southern tail in NGC 
4038/39. Since then, a lot effort has been spent on obser- 
vations of candidate objects for such tidal dwarf galaxies 
(TDGs). 

The distribution of gas and stars in tida l tails 
can be very di fferent l|Hibbard fc van Gorkoml 1 19961 : 
Hibbard fc Yunl [l999l: [H cithauscn & Walter] |2000| : 



Smith et all 1997: 



Hibbard et all |2000| ). Some tails 
show gaps which could be a consequence of local 
collapse or feedback associa ted with star formation 
jHibbard fc van Gorkoml Il996l ) and in some tails, there 
is an offset between t he stellar and the gaseous material 
llHibbard fc Yunlll999l ). Tidal tails can easily contain up to 
several 10 9 M^ of g a s and a similar mass i n stars (see e.g. 
Hibbard et all 1 19941 : IHibbard fc Yunl 1 19991 : IHibbard et all 



2 00l j) . The molecular complexes for ming in the tails can 



reach masses of a few 10 7 Mp, (see e.g. lHeithausen fc Walter! 
l200d : IWalter fc Heithau"senlll999D . 

The formation processes of TDGs and star clusters in 
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tidal tails seem to be linked, as tidal tails with TDG candi- 
dates seem to have few star clusters, while those with many 
star clusters have no TDG candidates (|Knierman et al.l 
2003). Both star cluster formation as well as TDG forma- 
tion are directly coupled to the local properties of the gas 
reservoir. The formation of TDGs might rely on ex ceeding 
a threshold in HI column density (|Duc et al.ll2000l ). much 
like the similar threshold for star formation in dwarf ir- 



regulars and low surface brightness galaxies 


Davies et al. 


1976 


Gallagher & Hunter 


1984l;|Skillmanlll987 


Taylor et al. 


1994 


van der Hulst et al.l 


19931; Ivan Zee et al] 


1997T). As the 



1986), TDGs should have a lower threshold than dwarf 
galaxies formed in isolation because TD Gs are made fro m 
recycled material of the progenitor discs (iDuc et al.l feOQO'), 
To the present day, t here are a r ound 60 well known 



TDG candidate ob j ects (ISchweizerl 1 19781; iMirabel et al 
19921 ; iHibbard et all 1 1994 lYoshida et al l Il994l; IDuc et al 



2000 



1997 ; IDuc fc Mirabell 119981 ; iHeithausen fc Walter 



Due et al 



2002: 



iMendes de Oliveira et al 
20011; IWeilbacher et all |2002|; iTempqrin et all bool 
Weilbacher et al.ll2003l;lMundell et al.ll2004l ; lBournaud et all 
2004l ; lAmram et al.ll2004r ). However, the quality and level of 



details of the observations still varies considerably for these 
objects, which also introduces some scatter in the reliability 
of their detections as TDG candidates. 

In general a TDG can be detected either while it is 
still embedded in the tidal tail or as a dwarf galaxy ac- 
companying an already merged system. In the first case, 
one will try to detect overdensities in gas and stars in- 
side the tidal tail, with the ultimate proof for a newly 
formed TDG being th at it kinematica lly decoupled from 
the surrounding tail (|Duc et all |2000T ) . However, corre- 
sponding signatures in velocity gradient o r dispersion can 
easily be masked by projection effects l|Bournaud et all 
120041 ') a nd so far only only o ne such object has been con- 
firmed (|Bournaud et al. ll2004~) . while for other observed ob- 
jects the velocity gradien ts are not large enough to as- 
sure self-gravity (see e.g. 



IMendes de Oliveira et all | 200l|; 
Hibbard et al.ll200ll;ISaviane et al.ll2004l ; lAmram et al.ll2004 



Bournaud et al.ll2004h . 



Alternatively, one can detect TDGs after their de- 
coupling from the tidal tail, i.e. to search the popu- 
lation of dwarf galaxies in merger systems for TDGs 

1 



Deeg et all 1 19981 ; iDelgado-Donate et all 12003(1 



IDuc fc Mirabel ( 19981 ) suggested to use a metallicity based 
selection criterion of «Zq/3 for TDGs. Given the higher 
metallicity of TDGs compared to usual dwarf galax- 
ies, TDGs would thus also deviate from the luminosity- 
metallicity relation of dwarf and giant galaxies (see e.g. 
Skillman et ai1ll989l ; iRicher fc McCalllll995l ; ISkillman etatl 



19971 ). 



Recently, a s ample of 13 TDG candidate objec ts has 
become available (|Weilbacher et al.ll2000l . |2002| , 120031 ). The 
majority of these TDG candidates has luminosities of 
Mb ~ —12 to —17 in B-band and a mean Ha luminos- 
ity of 2.2 xlO 39 erg s _1 . Their mean oxygen abundance is 
12 + \og(0/H) = 8.34 ± 0.2, consistent with material in the 
outer disc of spiral g alaxies like the Milky Way (see e.g. 
iNaab fc Ostrikerll2006M . All 13 objects seem to have internal 
velocity gradients, but higher resolution spectra would be 
required for confirmation. 



1.2 Theoretical Modelling of Tidal Tails 

The formation and evolution of tidal features in systems of 
interacting galaxies has been s tudied with theoretical models 
since the early 1 960s (see e.g. IPfleiderer fc SiedentopJ I l961 
iPfleidcrcr 1963]; iToomre fc Toomrel Il972l; lYabushitaf 11971 



IClutton-Brockl Il972al lbl; IWrightl Il972l ; lEneev et all Il973l ). 

Most of these models relied on restricted three-body sim- 
ulations. With the advent of more powerful computers, also 
more advanced algorithms for such particle simulations be- 
came available. Tree codes have been the method of choice 
for many merger studies, as they offer a wide dynamic 
range regarding their spat i al and mass resolut i on (see e.g . 
Appcl 1985; Jernigan 1985; Barnes fc Hutlll986l ; |Presslll986l; 
Jernigan fc Por ter 1989; BarnesHl990T ). Gas has been incor- 



porated into such tree code simulations using th e Smoothed 
Particle Hydrodynamics method, or short SPH ijLucvl 19771 ; 
iGingold fc Monaghanlll977l ; lBenj|l990l ; lMonaghanlll992r i. 

Self-consistent merger simulations have been used since 
the mid 1980s to study the formation of the remnant 



galaxy and its properti e s. Gen era l stud i es of merging sys- 
tems (see e.g. iBarnesI ll988bl R 1 19891 ; INaab et all Il999l ; 



INaab fc Burkertll2003l ) have been conducted as well as at- 
tempts to model specific observed merger systems (see e.g. 
Mihos et alll 19931; [Hibbard fc Miholl 19951; lMihos fc Bothunl 



19971 ; iMihos et alll 19981 : IDuc et al.ll2000l ; lBarnesll2004r i. The 



distribution of gas and stellar material in the tidal tails is 
of major importance to determine the encounter geometry 
of the observed system. However, those simulations repro- 
duced the tidal tails as such and were not particularly aimed 
at modelling the formation of bound structures in the tidal 
tails. 

Reproducing the star formation history of merging sys- 
tems has been another constraint for simulations. However, 
the star formation rate can easily be underestimated in re- 
gions where strong compression due to shocks in the en- 
counter occur, which also makes it difficult to model merger 
induced starbursts which are strong en ough to yield emis - 
sion comparable to observed ULIRGs |Mihos et all ll993T ). 
More refined star formation prescri ptions in the models are 
required to model such star bursts (Barnes 2004). 

Merger simulations have also been successful in un- 
derstanding the offsets between gas and stellar mate- 
rial observed in some merging systems llSmith et all 1 19971 : 
IHibbard fc Yunll 19991 ; THibbard et allboooj ). Rather than be- 
ing due to complex feedback phenomena in the tails, those 
offsets are due to dissipational effects during the interac- 
tion, which lead to a kinematical decoupling of the gaseous 
compon ent from the stellar one soon after the tidal tails are 
ejected (|Mihosl 120011 1 Hence such offsets are a kinematical 
consequence of some specific encounter geometries. 

The tidal debris of a merger can be used to probe 
the properties of the underlying dark matter halos 
llFaber fc Gallagherlll978l ; IWhiteil982l ; iNegroponte fc White! 
1 19831 ; I BarnesHl988bT l . Merger models in which the dark mat- 
ter halos are se t up as pseudo-i sothermal spheres (see e.g . 
lHernquisd[l990l) or NFW halos (|Navarro et al.lll996l . 119971 1 
can form significant tidal tails, provided that the ratio of es- 
cape velocity « e to circular velocity v c at ab out the solar ra- 
dius i n the progenitor galaxy is v e /v c ^ 2.5 l|Dubinski et all 
1 19991 ; ISpringel fc Whitelll999l ) . 

Most merger simulations concentrated either on the 
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properties of the remnant galaxies or on the global dis tribu- 
tion of mass in tidal tails. IBarnes fc Hernqu ist (1992) were 
the first to concentrate explicitly on the formation of struc- 
tures inside the tidal tails. They modeled one encounter, 
corresponding to one pair of tidal tails, in which they found 
23 clumps of material, possibly being TDGs or their pro- 
genitors. The resolution of these tails was limited by the 
low numbers of particles in the simulation, in total « 45000 
per galaxy with ~ 8000 of them for the gas component. The 
most massive clump of material in the tidal tails was re- 
solved with 350 particles and had a total mass of 4.18 x 10 s 
Mp; if the progenitor galax ies are scaled to the Milky Way 
ijBarnes fc Hernquist|[l992l ). 25% of the mass inside the ob- 
ject was in gas. The dark matter content of the object was 
< 5%. The distance from the remnant at which the most 
massive object forms varies between pure N-body simula- 
tions of the sa me encounter IIBarneslll992f) and simulations 
including gas (|Barnes fc Hernquistl Il992l . Il99rj ). While the 
stellar material which ends up in the TDG candidates gen- 
tly collapses once removed from the parent galaxy, the gas 
is driven into a thin 'ridge line' along the tail by convergent 
flows. The bound objects were thought to form in the stel- 
lar co mponent alone, then gas assemb les in these potential 
wells l|Barnes fc Hernquistlll992l . Il996l ). This formation sce- 
nario has been one of the two standard p aradi gms for the 
formation of TDGs. IBarnes fc Hernquistl (|l996h emphasize 
that their pure iV-body simulation showed similar objects 
than the one including gas and thus point out that "dis- 
sipative effects are not crucial in forming such structures", 
as the density fluctuations in the stellar component alone 
were sufficient to produce the TDG candidates. This should, 
however, be taken with care, as the source for the density 
fluctuations in the tail which then collapse and form TDG 
candidate objects are most prob ably particle noise and swin g 
amplification during the merger IBarnes fc Hernquistlll992f ). 
which means that at different resolutions and thus different 
particle noise levels, the fragmentation of a tidal tail possi- 
bly looks quite diffe r ent w hich hasn't been addressed since 
IBarnes fc Hernquistl (|l992T ). We will attempt to do so in sec- 
tion [3] 



lElmegreen et al.l (|l993l ) and lKaufman et all i|l994h sug- 
gested a different process for the formation of bound clumps 
in tidal arms. This scenario, in which gas plays the key role 
in the formation of TDGs is the other major paradigm es- 
tablished for the form atio n of TDGs. The simu lations of 
lElmegreen et alj (|l993h and lKaufman et al.1 (|l994h show the 
growth of large cloud complexes as massive as m 10 s Mq in 
the tidal arms. The contr ibution of gas to t he to tal mass of 
those objects was « 60%. lElmegreen et all (|l993h explained 
the formation of these objects with the high velocity dis- 
persion in the gas of interacting systems. A cloud forming 
from gravitational instabilities inside an ISM with increased 
velocity dispersion is likely to be more massive, as the local 
Jeans mass of the gas is increased as well. In addition, the 
dispersion intrinsic to the cloud should also be increased, 
allowing it to form stars with higher efficiency as the cloud 
is more stable against self-destruction. The conclusion of 



lElmegreen et al.l (I1993T ) and lKaufman et al.1 1| 19941 ) that gas 
is the key ingredient f or the formation of TDGs c learly con- 
tradicts the results of IBarnes fc Hernquistl l|l992T) . It should 
be noted that the numerical models used to arrive at these 



two results differ considerably in terms of the numerical 
methods as well as the initial conditions. 

iLi et al.l |2004l ) simulated the formation of globular 
clusters in tidal tails. Although they were able to reproduce 
the distribution and metallicity of globular clusters in ellip- 
ticals, the insight on TDG formation in those si mulations is 
rathe r limited because of the use of sink particles (|Bate et al.l 
1995) for modelling the collapse of gaseous structures into 
globular clusters. As every sink particle resembles one globu- 
lar cluster, no larger stellar dynamical systems like TDG are 
modelled. Stellar superclusters formed in mergers possibly 
merg e later and form ultracomp act dwarf galaxies (jKroupal 
1 19981 ; iFellhauer fc Kroupall2002l ) . 

Recently there has been some discussion about the con- 
ditions under which TDGs can be formed at the very tip of 
the tidal tail, a location in which they or their progenitors 
are frequently found in observed systems (note that not all 
TDG candida tes are located at the ti p of t heir ambient tida l 
tail, though). IBournaud et all |2003l ) and lDuc et ail (I2004T I 
have argued that there could exist two distinct formation 
mechanisms for TDG candidates: Local gravitational insta- 
bilities in the inner part of tidal tails lead to the formation 
of TDG candidates with masses of 10 7 - 1O 8 M while mas- 
sive (> 1O 9 M0) TDG candidates at the very tip of tidal 
tails form in a top down scenario in which the tidally ex- 
pelled gas accumulates at the tip of the tail and then under- 
goes gravitational collapse. Thi s process might be similar t o 
the process of co ll apsin g shee ts llBurkert fc Hartmannl 2004). 
IBournaud etall l|2003l ) and iDuc et al.1 l |2004h argue for a 
threshold in the radial extension of the progenitor galaxy's 
dark matter halo relative to the size of the progenitor's stel- 
lar disc. Only haloes with cutoff radii at least 10 times larger 
than the radial scale length of the disc should hence be able 
to form such massive TDG candidates at the tip of tidal 
tails. 

However, some care is advised in the interpretation of 
those simulations, as they have a rather coarse spatial reso- 
lution of the iV-body component of 5 kpc, which is why the 
gravitation of the gas component in IBournaud et al.1 (|2003h 
is only simulated in two dimensions, while the rest of the 
simulation is three dimensional. The 2D resolution is 150 pc 
at the cost of neglecting the extent of the gas in the tail 
perpendicular to the orbital plan e completely. This problem 
is overcome in IDuc et al.l l|2004l ). but still the resolution is 
only 390 pc in one simulation and even a factor of two worse 
for all other simulations. Both studies don't model the gas 
component with full hydrodynamics, but with the simpler 
sticky particles approach. 

In this paper we explore the origin of TDGs with high 
resolution simulations. In Section !^. ll we describe the initial 
conditions of our numerical models and in Section 12.21 the 
simulation code. The gravitational force softening used in 
our simulations is discussed in Section 12.31 The results of 
our simulations are presented in Section [3] where we discuss 
in detail the influence of numerical resolution on structure 
formation in the tidal tails in Section 13.11 The effects of 
gravitational softening on the results are discussed in Sec- 
tion 13.21 while the effects of gas in the tails are discussed in 
Section [3.41 In Section [3.51 we present detailed photometric 
and kinematical analysis of our most massive TDG candi- 
date. Finally, we summarize and conclude in Section [4] 
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2 NUMERICAL MODELS 



2.1 Initial Conditions 



Every galaxy in our simulations consists of an exponential 
stellar disc component, a stellar bulge and a dark matter 
halo. In addition, a subset of the simulations included also a 
gaseous disc. The galaxi e s have been set up with the method 
described by IHernquistl (|l993l ) and have been run in isola- 
tion to verify that they are in equilibrium for a time scale 
longer than the time until the first encounter in the merger. 
Note that for one model including a gaseous disk, we ver- 
ify the stability in detail in Section 13.41 We will only give 
a short summary of the parameters of our galaxy models 
and the system of units used . For d etails of the method, we 
refer the reader to IHernquistl l| 19931 ). We assume G = 1 and 
the galaxies are intrinsically scale-free. As a general rule, all 
quantities inside this paper are given in these intrinsically 
scale free units unless otherwise noted. The conversion to 
physical units is performed by scaling a galaxy to some sys- 
tem with known physical properties. If the progenitors are 
scaled to the Milky Way, the unit time is 13.1 Myr, unit 
length is 3.5 kpc, unit velocity 262 km s _1 and unit mass is 
5.6 x 10 10 M Q . 

The stellar disc has a mass of md = 1. The discs surface 
density profile is an exponential with scale length ra = 1 and 
vertical scale height zo = 0.2. The cut off radius of the disc 
is r max = 15rd and the vertical cut o ff is at z max = 10zo. 
The Toomre parameter (|Toomrell 19641 ) at the solar radius is 
Q = 1.5. A stellar bulge component with mass mi = 0.2m<i 
and scale length rb = 0.2rd is used. The cut off radius for the 
bulge component is r max ,b = 15r,j. The dark matter halo is 
realized as a pseudo isothermal sphere with mass mj, = 5.8. 
It has a core radius of r Ct h = 1 and a cut off radius of 
r C ut,h = 10. 

From the large sa mple of 112 merger simulations of 
iNaab fc Burkert] (|2003h with different mass ratios and en- 
counter parameters, we chose the encounter which produces 
the most pronounced tidal tail and hence is a most favor- 
ab le case for TDG forma tion. The initial conditions used 
bv lNaab fe Burkertl |2003l ) were consistent with initial con- 
ditions derived from cosm ological N-body simulations (see 
iKhochfar fe Burkerj200^ ). In the encounter we selected, the 
disc of the galaxy which produces the tidal tails is aligned 
with the orbital plane of the encounter, the galaxies have 
equal mass and the encounter is prograde. The galaxies are 
put on a parabolic orbit, the initial separation between their 
centers of mass is 30. The pericenter distance of the en- 
counter is r v = 2. Together with their masses, this deter- 
mines the orbit completely. 

In order to study the formation and subsequent evolu- 
tion of tidal tails in detail, the resolution of the progenitor 
galaxy from which the tidal tail emerges should be maxi- 
mized as far as possible in a simulation. The second galaxy, 
on the other hand, only acts as a perturber. It creates the 
tidal forces which strip the future tail material from the pro- 
genitor disc. The internal structure, dynamics and orienta- 
tion of the perturbing galaxy are irrelevant for the tidal tail 
as long as the tidal forces during the encounter are repro- 
duced accurately enough. In order to be able to maximize 
the resolution of the progenitor galaxy, we therefore only re- 
solve the perturbing galaxy comparatively coarse. This ap- 
proach allows the study of only one tidal arm per merger, 



but this in greater detail, which is the aim of this work. We 
have confirmed the validity of this approach in a test sim- 
ulation in which both galaxies were resolved at the same 
level. The matter which is found inside the tidal arm of one 
galaxy contains no particles of the perturbing galaxy. The 
tidal arm forms inside halo material of the galaxy to which 
also the stellar material of the tail once belonged. 

We set up the perturbing galaxy with a resolution of 
10000 particles for the stellar disc, 2000 particles for the 
bulge component and 20000 for the dark matter halo. Both 
the perturbing galaxy as well as the one which develops the 
tidal tail share the parameters for the mass and extent of 
the components as described above. However, the galaxy 
which develops the tidal tails is realized at much higher 
resolution than the perturbing one. The particle numbers 
for our different models are summarized in Table [T] The 
lowest resolution model, A, has already more particles in- 
side one galaxy than the e ntire simulation of the classical 
iBarnes fc Hernquist (1 19921) pa p er an d the follow up simu- 
lations of Barnes fc Hernquistl l|l996l ). As far as we know, 
our highest resolution simulation, model E, is the highest 
resolution simulation of a tidal tail published so far. 

2.2 Numerical Method 

The simulations were performed with our hybrid N-body / 
SPH code VINE. The code will be described in detail in 
a forthcoming paper (Wetzstein et al., in preparation), so 
here we will only give a short description of the numeri- 
cal algorithms used. VINE uses a binary, mutual nearest 
neighb or tree for the efficient c alculation of gravitational 
forces (|Appel|[l985l ; |Presslll986l ). A distant tree node j is 
accepted for the interaction list of particle i if their dis- 
tance Tij satisfies the condition r%j > hi + hj /6 where hi 
and hj are the physical size of the particle and node, re- 
spectively. For the tolerance parameter 8, we adopted a 
value of 0.7 for our simulations presented here. The time 
integration was carried out using a leapfrog integrator in 
combination with an individual time step scheme which al- 
lows the code to follow very high density contrasts in the 
simulated system very efficiently. The individual time step 
sche me is similar to the blocks tep scheme of the TREESPH 
code I Hernquist fc K atz 1989). The timestep i n +i of a par- 
ticle is computed usi ng the minimum of three different crite- 
ria: = min(ri y7i/|aj, T2/i/|v|, T2|v|/|a|) where h is the 
gravitational softening length, v the velocity and a the accel- 
eration of the particle. The Ti are tolerance parameters, for 
which we used T2 = 0.14 and T2 = 1. Our highest resolution 
simulations required time steps being a factor of more than 
32000 smaller than the maximum timestep in the scheme. 
The SPH s cheme in VINE follows the SPH formulation of 
Bcnz (1990j). The a rtificial v i scosity inclu ded the additional 
terms suggested by iBalsaral (|l990l . Il995i ) to reduce entropy 
generation in shear flows. Our choice of parameters of the 
code allowed us to keep the error in the total energy of the 
system lower than 0.3% in our JV-body simulations. 

VINE is fully parallelized with OpenMP. The gravita- 
tional force calculation, which is the most time consuming 
part of a simulation, scales linear up to 100 processors and 
the second most time consuming part, i.e. the SPH calcu- 
lation scales equally well. The simulations have been run in 
parallel on 32 processors of our SGI Altix supercomputer 
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Model 


-Ndisc 


■^bulge 


Nhalo 


N ga i 


Ntotal 


hstar 


hhalo 


A 


80k 


16k 


160k 


256k 


288k 


0.05 


0.071 


A2 


80k 


16k 


640k 


736k 


768k 


0.05 


0.035 


B 


160k 


32k 


320k 


512k 


544k 


0.039 


0.056 


C 


320k 


64k 


640k 


1024k 


1056k 


0.031 


0.044 


D 


640k 


128k 


1.280k 


2048k 


2080k 


0.025 


0.035 


D2 


640k 


128k 


160k 


928k 


960k 


0.025 


0.056 


E 


1280k 


256k 


2.560k 


4096k 


4128k 


0.019 


0.028 



Table 1. Particle numbers of the various components of the pure Af-body simulations. The 
columns are (from left to right): model name, number of disc particles, number of bulge particles, 
number of halo particles, number of particles inside the galaxy which develops the tidal tails, total 
number of particles inside the simulation, gravitational softening lengths of disc, bulge and halo, 
respectively. Numbers have been rounded to three digits behind the decimal point. 



at the University-Observatory, Munich, as well as the IBM 
Regatta supercomputer at the Rechenzentrum of the Max 
Planck Society, Garching. 

2.3 Gravitational Softening 

As we want to study the collapse of objects inside the 
tidal arms, the particle noise is an important factor 
( Barn es" fe Hernquistlll992h . An additional effect here is the 
gravitational force softening used. The optimal choice for 
the force softening, regarding its functional form as well as 
the particular softening scale used as a function of particle 
number have been a subject of intensive debate in the litera- 
ture (see e.g. iMerrittll 199q ; iRo mco 1997; lAthanassoula et all 
l200d : lDehnenl200i and references therein). We use the SPH 
spline kernel to soften the gravitational forces of all parti- 
cles. It has the advantage that outside the softening length 
it exactly resembles the Newtonian po tential. The be nefits 
of this choice ove r Plumm e r soft ening |Aarsethlll963l ) have 
been discussed by iDehnenl (1200 ll ). 

Regarding the choice of the softening lengths, there ex- 
ist no studies particularly dedicated to mergers of interact- 
ing galaxies. We have thus used empirical values for the soft- 
ening of the stellar disc in model A. For a merger with a 
disc of 80000 part icles, our choice for the softening length is 
h = 0.05 (see e.g. iNaab fe Burkertll2003() , which guarantees 
a stable thin disc after several rotation periods. In order to 
obtain softening lengths for the bulge and halo component 
of the galaxy, we have adopted the point of view of force 
softening being equivalent to giving a point mass a finite 
mass distribution. Particles with higher mass should thus 
have increased softening lengths as well. This idea was ap- 
plied to every component (and thus particle species with a 
given mass) of the galaxy by scaling the smoothing length of 
^soooo = 0.05 for our particles in the disc of model A (80000 
particles) with mass mgoooo = md/80000 to other masses 
according to h„ ew = /i80ooo(w-n eu ,/m8oooo) ( ~ 1/3) - The effects 
of two body relaxation and associated disc heating should 
hence be decreased, together with similar effects of dark 
matter halo particles penetrating the tidal tail at later stages 
of the simulation. The corresponding softening lengths used 



for the model galaxies are listed in Table [T] The soften- 
ing lengths of the components of the perturber galaxy have 
been adjusted in the same way, resulting in a fixed softening 
length of /idisc = 0.1 for the stellar disc, /ibuige = 0.1 for the 
bulge and /ihaio = 0.1797 for the dark matter halo. 

In any resolution study like the one presented here, soft- 
ening can affect the results. This can only be addressed by 
numerical experiments, i.e. actually varying the softening 
lengths and the particle numbers. We have carried out such 
tests in order to assure the validity of our results against 
softening effects, see Section \3. 2 1 for a detailed discussion. 



3 RESULTS OF THE RESOLUTION STUDY 

3.1 Influence of Increased Progenitor Resolution 

In order to study the effects of numerical resolution, we have 
conducted a sample of simulations with different resolutions 
of the same configuration, as described in Section [2TTJ Start- 
ing from a low resolution simulation with 80000 particles in 
the disc of the progenitor galaxy (model A), we increased the 
mass resolution up to a factor of 16 (m odel E). This should 
allow us to assess the hypothesis of iBarnes fe Hernquistl 
( 1992) that particle noise is one of the sources of the den- 
sity perturbations which lead to the gravitational collapse 
of objects. 

All the simulations have been run from t = until 
t = 200, at which the central part of the remnant galaxy 
has relaxed into a new equilibrium configuration. The first 
close encounter of the two galaxies occurs at t ~ 25 and the 
galaxies merge until t ~ 55. The time evolution of model 
A from t = 60 to t = 120 is shown in the upper panel of 
Figure [1] In this and the following particle scatter plots, we 
show only a randomly chosen subset of 40000 particles to al- 
low for better viewing. The tidal tail initially goes through 
a phase of rapid expansion. Later, its inner parts fall back 
onto the merger remnant, while the outer parts still expand 
and move further away from the remnant. 

Several fragments formed within the tidal tail and col- 
lapsed afterwards. They decoupled from the surrounding tail 
and became self-gravitating entities. For our analysis, we 




Figure 1. Time evolution of stellar material in model A (upper panel), and model C (lower panel), projected on the orbital plane. In 
model A, collapsed objects form in the tidal arm (last plot in upper panel), in contrast to model C (see last plot in lower panel). 



identify the objects with a friends-of-friends algorithm and 
check every object in the tail for self-gravity. We require the 
object to reach central stellar densities similar to those of 
observed dwarf galaxies, otherwise we don't consider the ob- 
ject as a possible candidate for a TDG. We adopt a stellar 
density threshold of po ^ 0.013 M^/pc 3 , corresponding to 
po,code ^ 0.01 in our code units. This choice is motivated 
by the low c entral density tail of observed dwarf galaxies 
(Mateo 1998), lowered slightly to yield a conservative lower 
limit. The formation of the objects in our model A is very 
simila r to those in the simulation of iBarnes fe Hernquistl 
l| 19921 ). We can confirm their result that dark matter does 
not significantly contribute to the overall mass of those ob- 
jects and find a fraction of ^ 3%. 

Model B evolves very similar to model A, hence we don't 
show its evolution in a separate figure. Several collapsed ob- 
jects form inside the tidal tails, but their distribution along 
the tidal tail is different from the one in model A. This in- 
dicates that results of simulations on TDG formation are 
affected by numerical resolution, as one would expect if par- 
ticle noise is affecting the gravitational collapse. 

Increasing the mass resolution by another factor of two 
seems to suppress the gravitational collapse, as can be seen 
from the lower panel of Figure Q] which shows the time evo- 
lution of model C. Apparently some regions in the tidal tail 



are more dense than others, but no collapsed objects as in 
the lower resolution models are formed. In order to study 
this result more quantitatively, Figure [2] shows the density 
inside the tidal arm as a function of distance from the center 
of the merger remnant for models A, B, C, D, and E. The 
upper panel shows the density at t = 120, the latest time 
shown in Figure [T] In addition, the lower panel shows the 
densities along the tidal tail at t — 200, a rather late phase 
of the tail evolution which allows us to assess the final stage 
of collapse along the tidal arm. We obtain the density in the 
tidal arm from the stellar components alone, as the contri- 
bution of dark matter to the total mass of the tidal arm is 
negligible. For every star particle, we compute its distance 
from the center of the merger remnant. The density of the 
particle is computed in a way similar to the SPH density 
calculation: For every particle, we search for the 50 nearest 
neighbouring particles and compute a mean density of the 
sphere containing those particles. For the plots in Figure [5] 
we create bins of equal mass in order to allow for a rea- 
sonable comparison between the different resolutions. Every 
point in the profiles of Figure [2] represents the maximum 
stellar density of the corresponding bin. The fixed mass in- 
side each bin is chosen as 1.25 x 10 -4 , corresponding to 10 
particles in model A and 160 in model E, respectively. 

Figure [2] shows that with increased resolution, no col- 
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Figure 2. Local stellar density pi oca i in the tidal tail vs. distance 
r from the center of the remnant at t=120 and t=200 for models 
A, B, C, D, and E. The shaded region at the top of both panels 
indicate a lower limit on central stellar densities of dwarf galaxies, 
see text for more details. 



lapsed objects form inside the tail anymore. Note that the 
innermost objects in the tail at t — 120 have fallen back 
into the remnant galaxy at t — 200 and hence are no longer 
present in the lower panel of Figure[2] The tidal tail in model 
B, which also forms collapsed objects, exhibits a similar evo- 
lution as model A. However, the growth of overdensities is 
delayed with respect to model A. 

As was already indicated in Figure [TJ the tidal tail of 
model C does produce overdense regions, but these hardly 
develop into collapsed objects. Following the trend of model 
B, the growth of overdensities is delayed further. Only at 
very late stages, some of these overdensities developed a high 
density contrast relative to their surroundings. 

Models D and E follow this trend as well. The mass as- 
sembly of overdense regions is delayed even further, and the 
associated regions do not gravitationally collapse to form 
distinct entities, such as the objects in model A or B. Com- 
paring the density evolution of the different models further, 
there are remarkable differences between models A to D. 
However, the density of the tidal tails in model D and E do 
not differ much and agree considerably well especially for the 
early stages of the tail evolution. We interpret this as a sign 
that for the time scale until t = 200 studied here, the sim- 
ulations start to converge somewhere around the resolution 
of model E. 

In order to assure that our result is independent of the 
simulation code used, we re-simulated models A and D, us- 
ing a different rando m seed for the setup and the Gadget2 
code (jSprinerell lioblh for the time evolution. We used an 
integration accuracy parameter of 0.02 and a tolerance pa- 
rameter for the tree accuracy of 10 -4 . The results of the 
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Figure 3. Same as Figure |21 for models A and D as well as the 
additional models AG and DG, which were set up with a different 
random seed and simulated with Gadget2. 



Gadget2 run are shown in Figure [3] Again we see that the 
density in the peaks decreases when resolution is increased. 

We thus conclude from this analysis that the gravita- 
tional collapse of objects inside a tidal arm, which has been 
reported fo r pure A-body simulations |Barnes fc Hernquistl 
1 19921 . 1 19961 ) , is strongly suppressed if the simulations have 
high enough resolution. The TDG candidate objects formed 
in models A and B are artifacts of the poor resolution in 
those models. Furthermore, only these two models are able 
to form TDG candidate objects with cen tral densities sim- 
ilar to those of observed dwarf galaxies (|Matealll998h . All 
other models shown in Figure [5] fail to reach such central 
densities in objects in the tidal arm. As we described in Sec- 
tion 12.11 the particular encounter which we study here is 
a most favorable case for the formation of extended tidal 
arms. Hence it is unlikely that other encounter parameters 
for the merger would work better to produce TDGs. We can 
conclude that pure A-body simulations of merging galaxies 
fail to generate collapsed objects in the tidal tails which re- 
semble dwarf galaxies, provided that the resolution of the 
simulations is high enough. 

3.2 Effects of Gravitational Softening on Tidal 
Tail 

In resolution studies like the one presented here the total 
number of particles is not the only parameter affecting the 
local gravitational potential. The other key parameter in 
this context is the gravitational softening length. Both, the 
softening length and the number of particles can cause two 
body relaxation effects such as artificial heating of an N- 
body system. For example, a system realized with few par- 
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Model 


hstar 


hhalo 


C 


0.031 


0.044 


CL4 


0.0078 


0.011 


CH4 


0.124 


0.176 


D 


0.025 


0.035 


DL4 


0.0063 


0.0088 


DL16 


0.0016 


0.0022 



Table 2. Gravitational softening lengths of the star and dark 
matter particles in simulations used to study softening effects. 
Apart from different softening lengths, models CL4 and CH4 arc 
identical to model C, DL4 and DL16 are otherwise identical to 
model D. In CL4, the softening lengths have been decreased by 
a factor of 4 compared to model C, while in CH4 they have been 
increased correspondingly. Model DL4 has a factor of 4 lower 
softening than D, DL16 has a factor of 16 lower softening than 
D. 



tides can have a very noisy gravitational potential if the 
softening lengths are small, or a rather smooth potential if 
the softening length is increased considerably. Especially if 
the gravitational collapse of small structures is to be stud- 
ied, both softening length and number of particles should 
ideally be explored to evaluate their effects on the results 
obtained from such simulations. As the exploration of this 
two parameter space is usually very expensive in terms of 
cpu time, such studies are rather rare. The usual approach 
is to resolve the system at hand with the highest particle 
number affordable under given cpu time constraints and to 
simulate this setup with a hopefully well chosen gravita- 
tional softening length. The latter is extremely dependent 
on the simulator's personal experience, as systematic stud- 
ies of force softening as a function of particle number are 
rare and usually only exist for certain well behaved mass 
distributions which in turn are often very different from the 

[Merritt] Il99rj; iRomeol 1 19971 ; 



system to be studied (see e.g 



lTheislll998l ; lAthanassoula et al.l 



20001 ; iDehn en 2001, and ref- 



erences therein). 

Nevertheless, we have explored the effects of different 
gravitational softening lengths for some of our models. The 
set of simulations used for this purpose is summarized in 
Tabled 

We have repeated the simulation of model C, but with 
a factor of 4 lower (model CL4) and a factor of 4 higher 
(model CH4) softening length. A similar set of models has 
been created for case D, with factors of 4 and 16 lower soft- 
ening lengths in models DL4 and DL16, respectively. Apart 
from the different softening lengths, all other parameters of 
those models have been inherited from models C and D, re- 
spectively. Figure [4] shows the density in the tidal tails of 
models C, CH4, CL4, D, DL4 and DL16 as a function of 
distance from the center of the merger remnant. 

Although different softening lengths lead to the forma- 
tion of different high density regions in the tidal arm, none 
of these regions reach densities of dwarf galaxies. We can 
therefore conclude that our particular choice of softening 
lengths in models A to E has not artificially suppressed the 




120 



Figure 4. Same as Figure [2] for models CH4, C, CL4, D, DL4 
and DL16 at time t = 200. 



formation of bound objects with central densities similar to 
dwarf galaxies. 



3.3 Influence of Progenitor Dark Matter Halo 

From the analysis presented so far, it is not clear which 
component of the progenitor galaxy, i.e. disc, bulge or halo, 
plays the key role in triggering the collapse of objects in the 
tidal arm at low resolutions. The bulge component is very 
unlikely to play an important role here, as the vast majority 
of the bulge particles is located inside the merger remnant 
rather than the tidal tail. So we need to address the question 
if the resolution of the stellar disc of the progenitor galaxy 
is important or if the resolution of the dark matter halo 
of the progenitor galaxy plays a crucial role in this respect. 
The halo component could cause perturbations in the stellar 
material of the tidal arm which then could grow and finally 
collapse. However, it could as well be the resolution of the 
stellar component itself which causes the collapse. 

In order to study the effects of both, we have created two 
additional models, A2 and D2. They are identical in setup 
to model A and D, respectively, but their dark matter halos 
have been realized with different particle numbers (see Table 
[T]). In model A2, the halo bas been sampled with 640000 
particles, while in model D2 the halo has been sampled with 
only 160000. Model A2 shares the resolution of the stellar 
disc with model A, but the dark halo is resolved with four 
times more particles than in A, thereby giving A2 the halo 
resolution of model C. Model D2 also shares the resolution 
of the stellar disc with model D, but its dark matter halo 
is resolved with eight times less particles, giving it the halo 
resolution of model A. 

Apart from the different disc / halo resolutions dis- 
cussed above, the simulations of models A2 and D2 shared 
the identical set of parameters as used in the other models. 
The density profiles of the tidal tails for models A2 and D2 
are compared with the results of models A and D in Figure 

El 

In the tidal tail of model A2, the density evolution of 
the outer tail is very similar to that of model A. In the inner 
region of the tail, however, instead of two overdensities in 
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Figure 5. Same as Figure l2l for models A, A2, D and D2 at 
t = 200. 



model A there forms only one overdensity in model A2. The 
peak densities are almost identical. The overall behavior is 
similar in the sense that dense objects are actually formed 
in both A and A2. 

Model D2 also evolves rather similar to D. There are 
only small differences between the two simulations. In model 
D2 there forms a density peak at r ~ 91 which is not 
present in D. Apparently, the formation of this overdensity 
was caused by the lower resolution dark matter halo of D2. 
But apart from this one region, the results of D and D2 are 
similar. Collapse of objects along the tidal tail is still highly 
suppressed. 

If shot noise from the dark halo plays an important role 
regarding the collapse of objects in simulated tidal tails, this 
should be visible in model A2 as well as D2. The tidal tail 
in D2 results from a stellar disc which is as well resolved 
as in model D, which did not form collapsing objects in its 
tidal tail. Thus if in model D2 collapsed objects would have 
formed in the tail, this would indicate that shot noise of the 
dark halo plays an important role. In addition, as model A2 
is otherwise identical to model A, its halo has been given 
the resolution of model C. So if the resolution of the halo 
would dominate the collapse of objects in the tidal tail, then 
this collapse should be effectively suppressed in model A2. 
However, model A2 still shows collapsed objects as in model 
A while at the same time the lower resolution halo of model 
D2 only triggers the collapse of one region in the tidal tail. 
The rest of the tail is very similar to model D. 

In summary, the high resolution dark matter halo in 
model A2 did not at all suppress the collapse of structures 
inside the tidal tail. The evolution was very similar to the 
model with low resolution dark halo, A. On the other hand, 
the high resolution stellar disc of model D2 did not exhibit 
collapsing objects like in model A, for example, although the 
disc in D2 resided in a low resolution dark matter halo. The 
evolution of the tidal tail was comparable to that of model 
D, apart from one additional overdense region which formed 
in D2. 

If we would have chosen a single softening length for 
all particles in the simulation, as is often done when using 
Plummer softening, this softening length would need to be 
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Figure 6. Surface density distribution of gas and stars in the 
progenitor galaxies for models Gl, G2, EG1 and EG2 (upper 
panel) and ratio of stellar to gas surface density in those models 
(lower panel). The region with r 30 kpc is affected by low 
particle numbers in the corresponding radial bins. Note that for 
ease of comparison, the model has been scaled to the Milky Way 
and quantities are given in physical units. 
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Figure 7. Initial radial profile of the Q parameter for the gas 
disk. Note that for ease of comparison, the model has been scaled 
to the Milky Way and quantities are given in physical units. 



large enough to keep the most massive particles, i.e. those in 
the halo, collisionless. This choice would probably result in 
larger softening lengths than those which we used, possibly 
suppressing local collapse to some extent. 

We conclude that the major effect leading to artificial 
numerical collapse of structures inside tidal tails of pure 
A^-body simulations is the resolution of the stellar disc of 
the progenitor galaxy which forms the tidal tail. With high 
enough resolution, collapse of such structures is strongly 
suppressed. The resolution of the dark halo has some in- 
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Model Ngas Mgas r gas /r st ar 



10° 
10"' 



Gl 


32k 


0.1 


1.0 


G2 


45k 


0.3 


1.0 


EG1 


32k 


0.1 


3.0 


EG2 


45k 


0.3 


3.0 



Table 3. Parameters of the gas simulations. The first column 
gives the number of SPH particles used, the second the total mass 
of the gas component and the third the scale length of the gas 
disc in units of the stellar scale length. 
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Figure 8. Circular velocity (upper panel) and surface density 
(lower panel) for the high resolution galaxy in model EG2 when 
evolved in isolation. Note that for ease of comparison, the model 
has been scaled to the Milky Way and quantities are given in 
physical units. The solid lines show the initial state, the dotted 
lines are at t = 327.5 Myr and the dashed lines at t = 589.5 Myr. 



fiuence on the formation of overdense regions, but the effect 
is small compared to that of the stellar disc resolution. 

3.4 Gas Dynamical Effects on the Formation of 
TDGs 

The resolution study of pure JV-body simulations presented 
in Section 13.11 resulted in the interesting fact that in such 
simulations, tidal dwarf galaxies are suppressed if the reso- 
lution of the simulations is high enough. 

As gas is observed in practically any tidal tail (see Sec- 
tion [T7T|, it could possibly be the lackin g ingredient in simu- 
lation s of tidal dwarf galaxy formation. iBarnes fc Hernquistl 
(1996) argue that gas would not play an important role in 
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Figure 9. Same as Figure[2]for models C, Gl, G2, EG1 and EG2 
at t = 120. 



triggering the collapse of structures inside a tidal tail. How- 
ever, their argumentation was based on their result that pure 
JV-body simulations of tidal tails, without any gas, all by 
themselves already lead to the formation of collapsed ob- 
jects and hence gas would not be required for this process. 
As we have shown that for our models, such collapse in pure 
iV-body simulations is artificial and induced by too low res- 
olution, it seems viable that we explore in how far gas could 
actually be important for the collapse of structures in tidal 
tails. 

In order to study the effects of gas on tidal tails, we take 
the setup of model C (see Table [1] for the details) and add 
a gaseous disc component using the SPH approach. Model 
C was the lowest resolution simulation in our sample which 
did not show collapsed objects like those in models A and 
B in the tail. In addition, including hydrodynamic effects 
in a merger simulation significantly increases the amount of 
CPU time required for completion. So model C is a good 
starting point for the study of gas dynamical effects when 
taking the increased computational expense into account. 

We have run a set of four models including gas, which 
is distributed in an exponential disc. Our models comprise 
a combination of discs which have the same scale length as 
the stellar disc or a scale length which is a factor of 3 larger. 
The mass of the gaseous discs is either 10% or 30% of the 
total disc mass, which is fixed at = 1 as before in model 
C. The scale height of the gaseous disk is zo, g = 0.02. The 
resulting combinations are listed in Table [3] The less massive 
discs in models Gl and EG1 have been resolved with 32000 
SPH particles, while the more massive ones in models G2 
and EG2 are represented with 45000 SPH particles. For the 
setup with the massive, extended gas disc in model EG2, 
we have matched the surfac e density of t he gas at the solar 
radius to the observations of lDamel (| 19931 ). adopting a value 
of E r0 = 10.2 Mq pc -2 . Scaling our galaxy to the Milky 
Way, we fix the corresponding surface density in our gas disc 
accordingly. Model EG2 hence is more likely to resemble the 
actual gas distribution in late type galaxies. Nevertheless, it 
is important to study the possible effects of different gas 
distributions in the progenitor disc on the distribution of 
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Figure 10. Time evolution of stellar component in model 61, projected on the orbital plane. No collapsed objects form. 



gas in the tidal arm. The initial surface densities of the gas 
components in our models are shown in Figure [6] 

In all models, the gas has been evolved with an isother- 
mal equation of state at temperature T = 12000i\\ The 
functional form of the gravitational softening of the SPH 
particles is the same spline kernel used for the iV-body par- 
ticles. The gravitational softening length of the SPH parti- 
cles is always identical to their SPH smoothing length. This 
approach prevents effects of artificial numerical c ollapse or 
stabilization against such (|Bate fc Burkertlll997l ). In addi- 
tion, we use the gravitational softening length of the stellar 
disc particles as a lower limit on the softening length of the 



SPH particles. Note that this can only lead to underestimat- 
ing the degree of collapse and not to artificially enhancing it. 
The remaining parameters of the simulation are as described 
in Section |2~21 

In order to demonstrate the stability of the models, we 
evolve the high resolution galaxy in model EG2 in isolation. 
In Figure [7] we show th e initial radia l profile of the Toomre 
stability parameter Q (|Toomrelll964l ) for the gas, modified 
for a composite disc made ou t of gas and stars ( Wang fc Silkl 
11994 iNaab fc OstrikerlboOfj ): 
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Here, k is the epicyclic frequency, a is velocity dispersion 
and E the surface density, indices g and * refer to gas and 
stars, respectively. Between 6 and 10 kpc, the disc is only 
marginally stable, but as we impose a lower limit on the 
gravitational softening of the gas particles, this slightly sta- 
bilizes the disc. We observe no local collapse during the evo- 
lution of the disc. The rotation curves and surface density 
profiles are presented in Figure [5] We show the profiles for 
t = 25 = 2.14T 1/2 and t = 45 = 3.85T 1/2 where T 1/2 is the 
rotation period at the half-mass radius. When the galaxy 
is scaled to the Milky Way, this corresponds to t = 327.5 
Myr and t = 589.5 Myr, respectively. The galaxy remains 



globally stable for these time scales. Note that the first en- 
counter of the two galaxies takes place at around t « 25, 
so by this time the perturbation due to the tidal interaction 
will be much stronger than any effects internal to the galaxy. 

In Figure [9] we show the stellar densities along the tidal 
tails in all four gas models. Model C is shown as well to allow 
for an easy comparison. The densities plotted in Figure[5]are 
obtained in the same way as those plotted in Figures [2][5] 
Note that only stellar material and no gas was taken into ac- 
count. Only models G2 and EG2 have developed collapsed 
structures in the tidal tail at t = 120. The central stellar 
densities inside those collapsed objects is higher than in the 
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Figure 12. Time evolution of stellar component in model EG2, projected on the orbital plane. Several collapsed objects form. 



low resolution pure JV-body models A and B. Those failed 
to reach densities similar to observed dwarf galaxies (shaded 
region in Figure [9}, while the objects in G2 and EG2 easily 
reach such central densities. To illustrate the formation se- 
quence of these objects, we show the evolution of the tidal 
tail in model Gl in Figures [TOllTTI and that of model EG2 
in Figures 1121131 The stellar material is shown in the first 
figures, i.e.[lO]for Gl and !12l for EG2, while the evolution of 
the gaseous component is shown in Figure [11] for Gl and 1 131 
for EG2. 

In model Gl, the gas is ejected from the progenitor 
disc together with the surrounding stellar material. It is 



then compressed into a thin line along the middle of the 
tidal tail. A significant fraction of the gas in the tidal tail 
flows back onto the remnant galaxy. Those regions in the 
tail which have initially high densities (see panel for t=50 
in Figure Hip are dispersed by the expanding tail so that at 
late stages (t=120) no density enhancements have survived 
in the tidal tail. However, in model EG2, the more extended 
setup of the gas disc leads to much more gas which is tidally 
ejected from the progenitor disc. Later in the evolution of 
the tidal tail in model EG2, the gas forms a transient bifur- 
cated structure (panels at t = 50, ...90 in Figure [T^J). The 
gaseous tidal tail extends about Ar = 20 further than the 
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Figure 13. Time evolution of gas component in model EG2, projected on the orbital plane. The logarithmic gas densities are color 
coded. Several collapsed objects form. 



stellar tidal tail. Both these features in the g as distribution 
are a l so present in observed tidal tails (see e. g. Hibbar d et al.l 
1 19941 : iHibbard fc Yunlll999l : Iffibbard et aHl2000L l200lf ). At 
later stages the tidal tail in model EG2 forms collapsing ob- 
jects (t ^ 80) in the mid region of the tail, in contrast to the 
pure iV-body model C and model Gl. 

Model EG2 forms at least 5 collapsed objects in the 
tidal tail, while in model Gl there exists none. Looking at 
the distribution of the objects in EG2 more closely, they are 
slightly offset from the center line of the tidal tail towards 
its leading outer edge. This is a result of the offset between 
stellar and gaseous material in the tail, which develops as 



a consequence of the different initial mass distribut ions and 
the dissipative nature of the gas (see [Mihos 20Cll], for de- 
tails). Apparently, the gas triggers local gravitational col- 
lapse in the tidal tail where its own density is maximal, as 
one would expect for gravitational instabilities in a gas dis- 
tribution. For the overall tidal tail, this leads to the collapse 
of both the stellar and gaseous components slightly offset 
from the center line of the stellar tail. After the gas has 
triggered the collapse, it settles in the centers of the objects 
formed. As our model does not include a numerical recipe 
for star formation and associated feedback, the gas forms 
very dense knots in the center of each stellar body. 
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Note that not only models G2 and EG2 formed col- 
lapsed objects. Also in model EG1, there forms a massive 
collapsed object. However, the object is located in the inner 
part of the tidal tail and the object falls back into the rem- 
nant galaxy where it is tidally disrupted. In the late phase of 
the tidal tail, as shown in Figure [§] the object has long dis- 
appeared. So it would be wrong to believe that massive gas 
discs in the progenitor galaxy, like in models G2 and EG2, 
are required to form TDGs. Also a less massive gas disc can 
lead to their formation, provided that the disc is extended 
enough. Also, the mere presence of a dissipational compo- 
nent does not lead to the formation of collapsed objects, as 
clearly indicated by model Gl. 

In order to quantify how much gas is torn out of the 
progenitor disc, we plot radial profiles of the cumulative gas 
mass for all models, see Figure [141 The profiles of the initial 
conditions are shown as well as the profiles at t = 120, the 
same evolved state of the tidal arm which is also shown in 
Figure [9] In model Gl, a mass of m g « 0.09 in gas is con- 
tained within the remnant and its immediate surroundings 
(r < 5). This corresponds to 90% of the total gas mass in 
this model. The distribution of the gas in the tidal tail is 
very smooth, no gaseous substructures are present in Gl. A 
similar smooth distribution can be seen in model EG1, al- 
though more mass is deposited further out in the tidal arm 
due to the more extended initial distribution. In model G2, 
the mass of the gas inside and close to the remnant is ~ 0.25 
while the remaining gas with a mass of w 0.05 is located out- 
side. Assuming an inner cut for the tidal tail of r — 20, the 
gas mass inside the tidal tail is m 9ita « « 0.03, more than a 
factor of 3 more than in model Gl. This clearly shows the 
effect of the gas distribution in the progenitor disc, as both 
Gl and EG1 have the same total mass and only differ in the 
corresponding distribution. In model EG2, much more gas is 
torn out of the progenitor disc than in any other of our gas 
models, as one would expect from the massive and extended 
disc in the progenitor. A mass of around 0.12 is located at 
r > 5, which corresponds to 40% of the total gas mass. As- 
suming again an inner cut for the tidal tail of r = 20, we find 
that a mass of ps 0.08 resides outside r = 20, corresponding 
to 27% of the total gas mass. Some of the density peaks in 
Figure[5]can also be found as jumps of the cumulative mass 
at the location of the object. 

W e conclude that in contradic tion to the previous find- 
ings of lBarnes fc HernquistJ (I1996T ). gas is the most impor- 
tant factor in the formation of collapsed objects inside the 
tidal tails of merging galaxies. Also the processes which lead 
to the formation of the offset between stellar and gaseous 
tails are in the end not responsible for th e suppression of 
TDG candidate objects, as suggested by iHibbard fc Yunl 
(1999), but in the contrary are closely linked to their forma- 
tion. The gas disc in the progenitor galaxy needs to supply 
enough gas into the tidal arm to trigger the onset of collapse 
there. The mere presence of a dissipative component is not 
enough. 
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Figure 14. Cumulative gas mass for models Gl, G2, EG1 and 
EG2 at t=0 (solid lines) and at t=120 (dashed lines) 
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Figure 15. Observable major axis rotation velocity (upper panel) 
and velocity dispersion (mid panel) profiles of the most massive 
TDG candidate, plotting gas and stars. The lower panel shows 
the ratio of rotation to dispersion. For details see text. 



3.5 Properties of a TDG Candidate 

The most massive TDG candidate object formed in our 
model EG2, located at (20,-48) in the t = 120 panel of Fig- 
ure d! has a mass of 0.062, corresponding to w 3.5 x 10 8 Mq 
if the progenitor galaxy is scaled to the Milky Way. About 
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Figure 16. Surface brightness profile of the most massive TDG 
candidate. The shaded region corresponds to twice the gravi- 
tational softening radius. Filled circles represent gas plus stars, 
while empty circles are for the old stellar population only. 



70% of its mass is in gas while there is no dark matter in- 
side the object. The velocity dispersion of the dark matter 
halo is too high to allow for their efficient capture in the 
substructures of the tidal arm. 

The object described above is our best candidate object 
for a TDG. However, if this system resembles observed sub- 
structures in tidal tails requires a more detailed study which 
infers directly observable quantities. In order to assess the 
properties of our candidate object, we convert the particle 
distribution into an artificial image which is then analyzed 
in a similar way as observations to avoid systematic differ- 
ences. For the details of the procedure, see iNaab fc Burker tl 
<|2003h . 

We find an ellipticity of our object of e ~ 0.4. In Figure 
1151 we show profiles of the rotation velocity, the velocity 
dispersion and the ratio of the two along the major axis. We 
plot gas and stars together, which is equivalent to assuming 
that all gas was converted into stars after the TDG has 
formed. Although there is a signature of rotation as well as 
a peak in the velocity dispersion, the magnitude of these are 
rather low: A maximum rotation velocity of 0.05 corresponds 
to 13.1 km s _1 for a progenitor like the Milky Way, while 
the maximum dispersion of 0.095 corresponds to 24.9 km 
s _1 . Such values might be missed in observations, especially 
when the object is possibly seen in less favorable projections. 

We obtained a stellar surface density profile of our TDG 
candidate, see Figure 1161 following the method used by 
INaab &: Truiillol l|2006h . The line represents the best fitting 
Sersic model, corresponding to an index of n — 10.8 and an 
effective radius of r e = 0.2, which corresponds to 0.7 kpc 
if the progenitor galaxy is scaled to the Milky Way. The 
high index is caused by the concentration of the gas towards 
the center. For the profile without the newly formed stars 
(open circles in Figure flu]) . we obtain index of n f» 2 and an 
effective radius of r e = 0.4, corresponding to 1.4 kpc. 

The stellar properties of our TDG candidate are sim- 
ilar to observed properties of local group dwarf galaxies 



rtMateol 19981') as well as dwarf galaxies in groups or clus- 
ters l|jerien et al.ll2000l ; iBarazza et alll2003l ; fde Riicke et alj 
2005). But given the uncertainty in the evolution of the gas 
in the shallow potential well, a more realistic model for star 
formation might well affect the results. We expect that val- 
ues for n and r e with and without gas outline the most 
extreme cases. 



4 SUMMARY AND DISCUSSION 

In this study we investigated the for mation of tidal dwarf 
galaxi es. S ince the early simulatio ns of lBarnes fc Hernquistl 
(11992ft an d lElmegreen et al.ll|l993h . not much work has been 
done regarding the detailed study of TDG formation. A well 
accepted model for the formation of TDGs has, up to now, 
still to be established. We have studied the importance of 
resolution and force softening on the results of numerical 
simulations of tida l tails. Although particle no ise has already 
been suggested by I Barnes fc Hernquistl (| 1992ft as an impor- 
tant effect in this respect, no study prior to this work had 
investigated in detail the effects of different resolutions on 
TDG formation. Performing an extensive resolution study 
of pure iV-body simulations, ranging from a total of 288000 
particles to 4.128 x 10 6 particles, we find that with increas- 
ing resolution the collapse of objects inside the tidal tail is 
subsequently suppressed. At a resolution level of 640000 par- 
ticles in the progenitor stellar disc, the tidal tail only devel- 
ops overdense regions which do not collapse within 200 code 
time units, corresponding to At = 2.62 Gyr if the progenitor 
galaxy is scaled to the Milky Way. This result is robust at 
even higher resolution. The collapsed objects formed in the 
low resolution models cannot be identified with dwarf galax- 
ies, either, as they don't reach central densities comparable 
to those of dwarf galaxies. For the initial conditions pre- 
sented here, we have shown that low resolution disk mergers 
lead to collapsed objects in the tidal arms whereas for higher 
resolution this collapse is suppressed. 

In order to assess the effects of force softening on such 
simulations, we have run a second set of simulations in which 
the gravitational softening length has been varied. Using 
these models, we have shown that the lack of dense, col- 
lapsed objects in the tidal tails for high enough particle 
numbers is not affected by our particular choice of softening 
lengths. Furthermore we have demonstrated that the res- 
olution of the stellar disc of the progenitor galaxy seems 
to be most important to suppress the artificial collapse of 
regions in the tidal tail. The resolution of the correspond- 
ing dark matter halo only plays a minor role in prevent- 
ing artificially collapsing objects. We hence conclude that 
whether or not we form collapsed objects mainly depends 
on the resolution of the init ial disk. The previous findings of 
iBarnes fc Hernquistl l|l996ft might originate from their low 
nuemrical resolution. Based on this study, we suggest the 
use of at least half a million particles for modeling the pro- 
genitor disc and twice that number for the dark matter halo 
in future simulations of tidal dwarf galaxy formation. 

We have studied the effects of gas in tidal tails using four 
different gas distributions in the progenitor discs. The pres- 
ence of a dissipative component alone is not sufficient for the 
formation of TDG like objects. However, if the gas reservoir 
in the progenitor disc is massive and / or extended enough, 
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the tidal tail contains enough gas to trigger the collapse of 
objects. The TDG candidates in our gaseous models easily 
reach central densities comparable to dwarf galaxies, in con- 
trast to their counterparts in low resolution pure JV-body 
simulations. Our findings favor a dissipati on supported sce- 
nario for T DG formation, as sugg ested bv lElmegreen et al.l 
l| 19931 ) and lKaufman et all (|l994T ). 

We have analyzed the most massive TDG candidate 
object in greater detail. It has a mass of w 3.5 x 10 s Mq if 
the progenitor galaxy is scaled to the Milky Way. The ob- 
ject contains no dark matter and 70% of its mass is in gas. 
There exists a clear signature of rotation but the maximum 
rotation velocity is 13.1 km s" 1 . The velocity dispersion pro- 
file shows a clear peak in the center, corresponding to 24.9 
km s _1 . The surface density of our candidate object can be 
reasonably fit with a Sersic profile of index n = 10.8 and 
effective radius of r e — 0.7 kpc. If we analyze only the stel- 
lar component, we find n « 2 and r e = 1.4 kpc. Our values 
for n and r e outline extreme cases, in which either no or all 
gas is transformed into stars. A detailed model for star for- 
mation might result in values somewhere in between. Stellar 
feedback could also expel some fraction of the gas from the 
system. We will address these interesting questions in future 
simulations. Given the uncertainties in the treatment of gas 
physics and star formation, we find reasonable agreement 
with observations. 

Our results demonstrate the importance of the gas com- 
ponent and the need for high resolution simulations in order 
to study the formation of TDGs in mergers. In this context, 
models should ideally also take star formation and stellar 
feedback into account. Such simulations including additional 
physical processes are beyond the scope of this paper, which 
aimed at studying more fundamental effects, including nu- 
merical resolution and force softening. More detailed models 
which take star formation and feedback into account will be 
investigated in future work. 
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